Long-lived 2s state of muonic hydrogen: population and lifetime 
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Ab initio study of the density- dependent population and lifetime of the long-lived (fip)2s and 
the yield of (/J,p)i s atoms with kinetic energy 0.9 keV have been performed for the first time. The 
direct Coulomb 2s — > Is deexcitation is proved to be the dominant quenching mechanism of the 
2s state at kinetic energy below 2p threshold and explain the lifetime of the metastable 2s state 
and high-energy 0.9 keV component of (/j,p)is observed at low densities. The cross sections of 
the elastic, Stark and Coulomb deexcitation processes have been calculated in the close-coupling 
approach taking into account for the first time both the closed channels and the threshold effects 
due to vacuum polarization shifts of the ns states. The cross sections are used as the input data in 
the detailed study of the atomic cascade kinetics. The theoretical predictions are compared with the 
known experimental data at low densities. The 40% yield of the 0.9 keV(^tp)i s atoms is predicted 
for liquid hydrogen density. 



Introduction. Exotic hydrogen-like atoms are formed 
in excited states and further evolution of their initial dis- 
tributions in quantum numbers and kinetic energy is de- 
fined by the competition of the radiative and collisional- 
induced processes during the atomic cascade. Muonic 
hydrogen (fJ,~p) is of special interest among the exotic 
atoms due to its simplest structure and possibility to in- 
vestigate a number of problems both the exotic atom 
physics and bound state QED. It is well- known, that 
2s state plays a particular role in {pTp) atom due to 
(2s — 2p) Lamb shift, Et, — 202.0 meV, and has no ana- 
log in hadronic atoms (pionic, kaonic, etc.) in which 
strong interaction leads to nuclear absorption or annihi- 
lation from this state. 

As it is expected, the measuring of the Et, in the 
so-called \ip Lamb shift experiment [l| permits to de- 
termine the mean-square charge radius of the proton 
with a relative accuracy of 10~ 3 and results in a bet- 
ter test of bound-state QED. The success of this exper- 
iment depends crucially upon the population, e 2 ° ng (per 
fip atom), and lifetime, 7"2° ng , of the 2s state at kinetic 
energies below El. This fraction is called the long- 
lived or metastable fraction of the 2s state as the Stark 
2s — » 2p transitions are energetically forbidden and the 
rate of two-photon transition to the Is state is negli- 
gibly small as compared with the rate of muon decay, 



4.55 x 10 5 s . The delayed K a X-rays induced 



during the collisions (see Q and references therein) can 
also occur but have not been observed [|| until now. 

A high-energy component of (np)i s with kinetic energy 
~ 0.9 keV was recently discovered 0, [|| from the anal- 
ysis of the time-of-flight spectra (at low gas pressures 
Ph 2 = 16 and 64 hPa) and attributed to non-radiative 
quenching of the long-lived 2s state due to the formation 
of the muonic molecule in a resonant (np)i s + H 2 colli- 
sion and subsequent Coulomb deexcitation of the {jppp) + 
complex (see and references therein). However, a theo- 



retical estimation of the quenching rate in the framework 
of the side path model [f| gives ~5x 10 10 s _1 at liquid- 
hydrogen atom density (LHD=4.25 xlO 22 atoms/cm 3 ) 
that is about an order of magnitude smaller than the 
value 4.4t 2 ;g x 10 n s -1 extracted experimentally @. 

Contrary to this assumption, in our recent paper 
we suggested that observed collisional quenching of the 
metastable 2s state and high-energy (/xp)i s component 
can be explained by the direct Coulomb deexcitation 
(CD) process, 
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however, the problem of the (fJ.p)2s metastability has not 
been the subject of this paper. 

In the present study we confirm our previous predic- 
tions 7] by the ab initio quantum-mechanical calcula- 
tions of the cross sections for the elastic scattering, Stark 
transitions and CD in the framework of the improved 
close-coupling approach (CCA) taking for the first time 
into account both the closed channels and the threshold 
effects due to the vacuum polarization shifts of the ns 
states. The calculated cross sections of the above men- 
tioned processes have been used in the detailed calcula- 
tions of the atomic cascade kinetics. The population and 
lifetime of the metastable 2s state, and the yield of the 
hot (0.9 keV) (np)i s component at wide densities range 
(10" 9 - 1) LHD are predicted. 

Improved close-coupling approach. The close- 
coupling approach has been applied [8| by the authors 
to describe the elastic scattering and Stark transitions in 
the exotic atom- hydrogen molecule collisions. Recently 
0, H, OH HHi the unified treatment of the elastic scat- 
tering, Stark transitions and CD in the collisions of the 
excited exotic (muonic, pionic and antiprotonic hydro- 
gen) atoms with hydrogen ones have been performed in 
the framework of the CCA. Here, we briefly remind the 
main assumptions and outline of the approach. 

Since (fi~p) is neutral and in the low-lying states much 
smaller than target molecule, the distortion of the tar- 
get electron structure during collision can be neglected 
and moreover the exotic atom collision with the hydro- 
gen molecule can be approximately treated as the exotic 
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atom-hydrogen atom scattering. These assumptions are 
supported in particularly by a good agreement of the ex- 
perimental data [h| and our theoretical predictions [HI 
for the kinetic energy distribution of n~p atoms at the 
instant of nuclear absorption at LHD. Further, since vac- 
uum polarization shifts of the ns states are about two 
order of magnitude more than both the fine and hyper- 
fine splittings, the muonic atom states are described by 
the non-relativistic hydrogen-like wave functions with the 
energies of ns states shifted (due to vacuum polarization) 
relative to the degenerate nl states (1^0). 

Thus the total wave function of the system (/xp + H) is 
expanded in terms of the basis states constructed as ten- 
sor products of the free muonic and hydrogen atom wave 
functions and angular momentum wave function of their 
relative motion. In the present study we use the basis 
sets in which both the open and closed channels are in- 
cluded. The coupled second order differential equations 
are solved numerically by the modified Numerov method. 
The differential and integral cross sections for individual 
nl — » n'V transitions have been calculated for the princi- 
pal quantum number n — 2 — 8 and kinetic energies that 
are of interest in the detailed cascade calculations. 

According to our study, the closed channel effects are 
more pronounced for the processes in the lower states 
(n = 2 — 5) at very low collision energies especially near 
and below ns — np thresholds. Here we produce some of 
our results for the 21 — > 21' (elastic and Stark scattering) 
and CD 21 — > Is transitions (see figs.l, 2) where the 
effects of the closed channels are the most significant and 
extremely important for the problem of metastable 2s 
state under consideration. All presented results below 
2p threshold and also the rates for the CD 2s, 2p — > Is 
transitions above 2p threshold have been calculated for 
the first time. 

Figures 1 and 2 show the (np) n =2+ H collisional rates, 
Xni^n'i' = LHDa n i->n'i'V, of the elastic and Stark scat- 
tering and CD at LHD as a function of the laboratory 
kinetic energy. The calculations have been performed 
for three variants of the basis set: the minimal basis 
(n = 1,2), the extended basis (n < n max = 7) and 
(for comparison) the basis without closed channels at all. 
At kinetic energy below 2p-threshold only a few lowest 
partial waves contribute to the elastic 2s — > 2s and CD 
2s — > Is cross sections and, as it is seen from Fig. 1, the 
proper description of the collisional processes is impossi- 
ble in the subspace of the open channels. The inclusion 
of the nearest closed 2p state into the basis leads to the 
tremendous increase both the elastic 2s — > 2s and CD 
2s — > Is rates (about 10 3 and 10 6 times, respectively). 
The further extension of the basis permits to achieve the 
convergence when the closed states of muonic hydrogen 
up to n max = 7 were included. Above 2p threshold the 
closed channel effect is practically negligible for the elas- 
tic and Stark scattering and is less pronounced than be- 
low 2p threshold for CD (see Fig. 2). 

The results can be qualitatively explained by the dif- 
ferent nature of the elastic/Stark scattering and CD pro- 
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FIG. 1: Fig. 1. Collisional rates at LHD for the elastic 
2s — > 2s scattering (solid lines), Stark 2s — > 2p transition 
(dotted line) and CD 2s — > Is (dashed lines) vs. the labora- 
tory kinetic energy. The thin lines show the results obtained 
with the minimal basis (n = 1,2), the thick lines - with the 
extended basis, including all states of (fJ,~p) up to n max = 7. 
For comparison the dash-dotted line shows the rate of the CD 
2s — > Is (multiplied by 10 4 ) calculated without the closed 
channels. The vertical line indicates the 2s — 2p threshold 
energy (0.435 eV). 
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FIG. 2: Fig. 2. The same as in Fig. 1 but for the elastic 
2p —* 2p (solid line), Stark 2p —* 2s ( dashed line) and CD 
2p — > Is (dash-dotted thin and thick lines) processes. The 
energy is referred from 2p-threshold. The radiative 2p — Is 
rate is shown with horizontal line. 



cess. The elastic and Stark scattering is defined by the 
interaction in a range about a few atomic units (a.u.) and 
the high partial waves give the main contribution to the 
cross sections at energies above ~ 1 — 2 cV. So, the de- 
tails of the short-range interaction (at R < 0.1 a.u.) are 
less important for description of these processes exclud- 
ing the scattering at very low energies. In contrast, the 
CD process is accompanied by the large energy release 



3 



(tens and hundreds eV) and occurs at substantially more 
small distances (about a few units of the /ip Bohr radius, 
n 2 o |ll ) so the details of the short-range interaction and 
correspondingly the closed channel effects become more 
important. Besides, the number of partial waves involved 
in the CD process is always essentially smaller than in the 
elastic and Stark scattering because the large centrifugal 
barrier for higher partial waves prevents the exotic atom 
from reaching the interaction range relevant for the CD 
process. 

Since the ratios between the rates of the collisional pro- 
cesses are independent from the target density, the useful 
observations can be derived from the results presented in 
Figs. 1,2 without the atomic cascade calculations. The 
CD 2p — > Is is strongly suppressed in comparison with 
the other cascade processes and can be neglected in the 
kinetics calculations. In contrast, the rate of the CD 
2s — > Is at kinetic energy about a few tens meV becomes 
comparable with the rate of the elastic 2s — > 2s scatter- 
ing and therefore the CD 2s — > Is can occur quenching 
the metastable 2s state. It is clear, the formation of the 
metastablc 2s state can not lead to the fast CD 2s — > Is 
transition before the muonic hydrogen is thcrmalized. Fi- 
nally, the comparison of the CD 2s — ► Is rate with the 
one of muon decay permits to conclude that lifetime of 
the metastable 2s state at densities less than ~ 10~ 6 
LHD is mainly defined by the muon lifetime while at 
more higher densities the situation is quite different and 
the detailed kinetics calculations arc needed. 

Atomic cascade and the initial (n, I, E)- distribu- 
tions. The cascade in the exotic atoms (see [HI [HI and 
references therein) is conventionally divided into classi- 
cal (for n > 9) and quantum-mechanical (for n < 8) do- 
mains. In the classical domain the results of the classical- 
trajectory Monte Carlo calculations (for details see |14l |) 
of the elastic scattering, Stark mixing and CD processes 
have been included with the molecular structure of the 
target taken into account. The cross sections of the ex- 
ternal Auger effect have been calculated in the semiclas- 
sical approximation through the whole cascade. In the 
quantum-mechanical domain (at n < 8) the differential 
and integral cross sections for all nl — ► n'V transitions 
have been calculated in the present version of the CCA. 

The initial distributions in the quantum numbers (n, I) 
and kinetic energy E of the exotic atom at the instant of 
its formation are very important for the problem under 
consideration, especially at very low target densities rel- 
evant for the Lamb shift experiment. In the simplest 
picture of the exotic atom formation (usually used in 
the cascade calculations) the initial principal quantum 
number is fixed at no ~ 14 (for muonic hydrogen) and 
the statistical Z-distribution is assumed. More elaborate 
studies [l5|, [lfl [TtJ taking into account the molecular 
structure of the target result in the initial n-distribution 
with the sharp peak at lower values (for muonic hydrogen 
at no = 11) and non-statistical ^-distribution. 

Since in the limit of the lowest densities (< 10~ 7 LHD) 
the atomic cascade is mainly defined by the muon lifetime 



and the rates of the radiative transitions, the informa- 
tion about the initial E- and (n, ^-distributions can be 
obtained from the K A-ray yields and the kinetic energy 
distribution of muonic atoms in Is state. In the present 
study this knowledge was used to choose the initial dis- 
tribution parameters. We accept here the Gaussian dis- 
tribution centered at no = 11 as n-distribution and 
the modified statistical Z-distribution (at each n value): 
~ (21 + l)e~ ai ( 2 ' +1 ) with a one fitting parameter a/. As 
for the initial kinetic energy, we use here two-exponential 
distribution with the parameters fitted to reproduce the 
integrated kinetic energy distribution of (np)\ s extracted 
[|[ from the analysis of (up) diffusion time at gas pressure 
Ph 2 =0.0625 hPa. 
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FIG. 3: Fig. 3. The density dependence of the relative K a 
yield in muonic hydrogen for different variants of the initial 
/-distribution: statistical (dashed line) and modified statisti- 
cal (solid line for «i=0.08 and dash-dotted line for Q;=0.1) 
distributions. The experimental data are from [3J (asterisks), 
pl| (diamonds) and (open circles). 



The calculated relative K a A-ray yields are shown 
in Fig . 3 in comparison with the experimental data 
[S E3, EH ■ The theoretical results illustrate the effect 
of the initial ^-distribution at fixed n-distribution. It 
is seen that statistical ^-distribution results in contra- 
diction with the experimental data at densities < 10~ 6 
LHD. These differences can not be explained by the pos- 
sible uncertainties in the rates of the collisional processes. 
Contrary to this, the modified statistical distribution 
(a; =0.08) for the first time leads as a whole to excellent 
agreement between theory and experimental data [|| . At 
densities higher than ~ 10~ 4 LHD the collisional transi- 
tions and Stark mixing are more efficient and the initial 
Z-distribution is practically forgotten during the cascade. 

Metastable 2s fraction: population, lifetime and 
yield of hot Is component. The population e^™ 8 of 
the metastable 2s state and the yield e^* of the 0.9 keV 
(up) is atoms (due to CD 2s — * Is below 2p threshold) 
calculated in the present version of the atomic cascade 
are shown in Fig. 4 for two variants of the initial Z- 
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distribution: statistical (a; = 0) and modified statistical 
( ai = 0.08). 




FIG. 4: Fig. 4. Density dependence of the population e^T 6 
(solid lines) and the yield eij* (dashed lines) of 0.9 keV {jip)i s 
atoms obtained in ab initio cascade calculations with statis- 
tical (thin lines) and modified statistical (ai = 0.08) (thick 
lines) initial /-distributions. The corresponding experimental 
data are from [|J (filled triangles) and [j| (circles). 

The calculations with the modified /-distribution pre- 
dict higher populations e^" 8 as compared with the ones 
obtained with statistical /-distribution and also with the 
populations derived from the experimental data analy- 
sis (for details see d, [1]) at densities below 10~ 4 LHD. 
According to our study, the metastable fraction is about 
1 % below 10~ 7 and slowly increases in the density range 
10~ 7 - 10~ 2 LHD to - 4%. Above 0.01 LHD the e 1 ™ 8 
grows faster reaching ~ 40% at LHD. This behavior and 
growth of the metastable fraction above 0.01 LHD is a 
consequence of the dominant role of the elastic scattering 
21 — > 21 and Stark 2p — > 2s transition over the radiative 
2p —> Is transition (see Figs.l and 2). It is worth not- 
ing that density dependence of the e^s" 8 obtained in the 
present study with the statistical I- distribution are sim- 
ilar to the one obtained earlier by Jensen and Markushin 
[l4| only below 10~ 4 LHD. At higher densities our calcu- 
lations do not confirm the sharper increase of the popula- 
tion up to 65% at LHD predicted in [14|. The theoretical 
prediction for the e^" 8 is in an agreement with the exper- 
imental result [5j obtained in measurements of 0.9 keV 
(m p)is a t low gas pressures pn 2 = 16 and 64 hPa. 

The calculated for the first time yields 6^°* of the hot 
(0.9 keV) fraction (/xp)i s (see Fig. 4) at densities below 
10~ 4 LHD have a quite different density dependence as 
compared with the one of the e^™ 8 - In the density range 
10~ 7 — 10~ 4 LHD, according to our calculations, the hot 
fraction increases from 0.01% to ~ 2%. Above 10~ 4 LHD 
the values of the and e!, " 8 is practically coincide. 
Such a behavior can be explained by the growth of the 
thermalized fraction of the (/ip)2s atoms when the den- 
sity increasing. At very low densities (below 10~ 7 LHD) 



the thermalized fraction of the (np)2s below 2p threshold 
is negligibly small and the yield of the hot fraction is less 
than 0.01%. At higher densities it sharply increases with 
the density increasing and reaches ~ 1% at 10~ 5 LHD. 
Above ~ 10~ 4 LHD the (np)2s atoms (below 2p thresh- 
old) are slowed down to low energies and then quenched 
in collisions by the direct 2s — > Is CD and as it is demon- 
strated in Fig. 4 the population el,™ 8 and the yield of the 
hot fraction 6^°* are equal to each other. According to 
our study, we predict that the yield of the 0.9 keV (pp)i s 
atoms and the yield of the protons with kinetic energy 1 
keV at LHD is about - 40%. 




FIG. 5: Fig. 5. The density dependence of the lifetime, r 2 ° ns , 
calculated with the extended basis set with n max =7 (solid 
line) and with the minimal basis set including only closed 2p 
state (dashed line). The experimental data are from [jj. 

The density dependence of the lifetime tv,™ 8 calculated 
with the extended basis set with n max = 7 in comparison 
with the result obtained with the basis including only the 
closed 2p state is shown in Fig. 5. The extended basis cal- 
culations are in a good agreement with the experimental 
data [f|. According to our study, at densities below 10~ 7 
LHD the lifetime of the metastable 2s state is mainly 
determined by the muon lifetime. At the density range 
10~ 7 — 10~ 5 LHD we observe the competition between 
the muon decay and the non-radiative quenching of the 
metastable 2s state due to the direct CD 2s — > Is. Above 
10~ 5 — 10~ 4 LHD the quenching of the long-lived 2s state 
is entirely determined by the CD process. 

Conclusion. For the first time the cross sections 
for the (pp) + H scattering process - elastic, Stark 
and Coulomb deexcitation - are calculated in the close- 
coupling approach including the closed channels. These 
cross sections are used as the input data for detailed 
study of the atomic cascade kinetics in wide density range 
(10~ 9 — 1) LHD. The results of our treatment allow us 
to assert that the direct Coulomb 2s — > Is deexcitation 
(calculated for the first time both above and below 2s— 2p 
threshold) is the dominant quenching mechanism of the 
long-lived 2s fraction and explains the experimental value 
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of its lifetime. The calculated population and lifetime of 
the metastable 2s-state and also predicted yield of the 
hot component (fip)i s with kinetic energy ~ 0.9 keV can 
be very useful for the choice of the optimal conditions for 
the fip Lamb shift experiment and other experiments. In 
particular, it will be very interesting to obtain experi- 
mental data about the density dependence of the yield of 
the hot component (np)u with kinetic energy ~ 0.9 keV 



or the yield of the protons with kinetic energy 1 keV. 
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